Effective chiral restoration in the p'-meson in lattice QCD 
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In simulations with dynamical quarks it has been established that the ground state p in the 
infrared is a strong mixture of the two chiral representations (0,1) + (1,0) and (1/2,1/2);,. Its 
angular momentum content is approximately the '^ Si partial wave which is consistent with the 
quark model. Effective chiral restoration in an excited p-meson would require that in the infrared 
this meson couples predominantly to one of the two representations. The variational method allows 
one to study the mixing of interpolators with different chiral transformation properties in the non- 
perturbatively determined excited state at different resolution scales. We present results for the first 
excited state of the p-meson using simulations with n/ = 2 dynamical quarks. We point out, that in 
the infrared a leading contribution to p' — p(1450) comes from (1/2, 1/2)^, in contrast to the p. Its 
approximate chiral partner would be a /ii(1380) state. The p' wave function contains a significant 
contribution of the '^D\ wave which is not consistent with the quark model prediction. 

PACS numbers: 11. 15. Ha, 12.38.Gc, ll.30.Rd 



1. Introduction. A central question for QCD is how 
both, confinement and chiral symmetry breaking, are in- 
terrelated and influence the mass generation of hadrons. 
Chiral symmetry is dynamically broken in the QCD vac- 
uum and this phenomenon is the most important for the 
mass origin of the ground state hadrons. It is believed, 
that a coupling of valence quarks with the quark conden- 
sate of the vacuum is responsible for the large constituent 
mass of quarks at low momenta. This large mass makes 
the problem effectively non-relativistic and the ground 
state p is a '^ 5*1 state according to the quark model lan- 
guage [l|-y| • Traditionally the excitation spectrum is also 
described with the quark model and, e.g., the first excited 
state of the p-meson, p(1450), is believed to be the first 
radial excitation of the ground state, i.e., the ^Si state 
ii- 

The empirical spectrum of the highly excited hadrons 
in the light quark sector exhibits patterns of parity dou- 
blets. It has been suggested that these patterns signal ef- 
fective restoration of both SU{2)l x SU{2)r and U{1)a 
symmetries in excited hadrons (5|, |6| (for a review see 
[7|). This conjecture would imply that the mass genera- 
tion mechanism in highly excited hadrons is very different 
and the quark condensate of the vacuum is of little im- 
portance. It would also imply that the constituent quark 
model language is not adequate for highly excited states. 

Such an interpretation of the spectroscopic patterns is 
not unique and it is possible to imagine some alternative 
explanations of existing symmetry patterns [8l4ll . Ef- 
fective chiral and ^7(1)^ restorations have a strong pre- 
dictive power. Some of the chiral partners for existing 
highly excited states are missing and one obvious way 
to confirm or reject effective chiral and U{1)a restora- 
tions is to find experimentally missing chiral partners or 
establish their absence jl2|. This is a difficult experimen- 
tal task, though. Another empirical approach to answer 
this question is to find alternative experimental signa- 
tures that would correlate with the spectroscopic pat- 



terns. Indeed, the effective chiral restoration would re- 
quire that the states with approximate chiral symmetry 
must almost decouple from pions and their diagonal axial 
coupling constants must be small Il3!,|l4|. It is difficult, if 
not impossible, to measure experimentally diagonal axial 
coupling constants for highly excited hadrons. The effec- 
tive chiral restoration also requires that the states which 
are assumed to be in approximate chiral multiplets, have 
a small strong decay coupling constants into the ground 
states and a pion. The analysis of empirical decays of 
excited nucleons shows that indeed all those excited nu- 
cleons that are approximate parity doublets have a very 



small decay coupling constant [15|. There are no other 
experimental tools that would help us to clarify this im- 
portant question beyond reasonable doubts. To resolve 
the issue one needs direct information about the hadron 
structure, which can be supplied in ab initio lattice sim- 
ulations. 

A first attempt to address the problem on the lattice 
was Ref. [l6| . The low-lying states are dominated by the 
near zero modes of the Dirac operator, which are directly 
related to the quark condensate of the vacuum. The role 
of these modes at small Euclidean times of the two-point 
correlation function is insignificant. This part of the cor- 
relator is dominated by the highly excited hadrons. This 
observation is consistent with effective chiral restoration 
but not sufficiently clear. A conclusive evidence would 
require to extract the high-lying chiral partners and see 
what the role of the near zero modes for the splitting of 
these states is. This is a difficult task, however. The 
other way is to measure the axial coupling constants of 
excited states [17|. For that purpose one would need 
reliable plateaus in two- and three-point functions for 
excited states near the chiral limit, which is not easy. 

Here we suggest an alternative approach to the prob- 
lem. Using the variational method 18] and a set of 
interpolators that scan a complete set of chiral repre- 
sentations for a given meson, we can study couplings of 



different interpolators to a given meson at different res- 
olution scales [19J. This method has successfully been 
applied to the ground state p-meson and the analysis has 
revealed that the ground state in the infrared is a strong 
mixture of the two possible representations (0, 1) + (1,0) 
and (1/2, 1/2)6 of SU{2)l x SU(2)r. Chiral symmetry 
is strongly broken in the state and this state is approx- 
imately a ^Si partial wave which is in agreement with 
the quark model language. Here we use this method for 
the first excited state of the p- meson, p(1450). We are 
able to present for the first time a clear evidence that in 
the infrared the structure of this state is different from 
the quark model prediction (^Si): it belongs predomi- 
nantly to the (1/2, l/2)b representation, which indicates 
a smooth onset of effective chiral restoration. 

2. Elements of the formalism. There exist two differ- 
ent local operators with the p-meson quantum numbers 
/, J^'^' — 1,1^^, the vector current, q^^rq, and the pseu- 
dotensor "current" , qa^'^rq. It is well established, both in 
quenched and dynamical (see, e.g., 19l-l24|) lattice sim- 
ulations, that both the ground state p-meson as well as 
its first excited state can be created from the vacuum 
by either of these operators. These two operators have 
distinct chiral transformation properties with respect to 
SU{2)l X SU{2)r 0,0, III. The vector current belongs 
to the (0, 1) + (1, 0) representation and its chiral partner 
is the axial-vector current qj^j^Tq, which creates from 
the vacuum the axial- vector meson oi (/, J^'~' = 1, 1^^). 
The pseudotensor interpolator transforms as (1/2, 1/2)^ 
and its chiral partner is the operator e^^^qa^'^q that cre- 
ates from the vacuum the /ii-meson (/, J^'^ = 0, 1^ ). 

Assuming that chiral symmetry is not broken both ex- 
plicitly and spontaneously, there would be two different 
groups of p-mesons. The first group would belong to the 
(0, 1) -|- (1, 0) representation, would couple exclusively to 
the vector current (or to the non-local operator with the 
same chiral transformation properties) and each member 
would be mass degenerate with its corresponding axial- 
vector partner ai. The p- mesons of the second group 
would transform as (1/2, 1/2);,, could be created from 
the vacuum only by the pseudotensor operator (or by 
the non-local operator with the same chiral transforma- 
tion properties) and would be systematically degenerate 
with their corresponding partners hi . A total amount of 
p-mesons in the spectrum would coincide with the com- 
bined amount of ai- and ft, i -mesons. 

Chiral symmetry breaking in the vacuum implies that 
the p-mesons are mixtures of both the (0, 1) + (1, 0) and 
(1/2, l/2)b representations and can be created by both 
the vector and the pseudotensor interpolators. Effective 
chiral restoration in highly excited p-mesons would re- 
quire that some of them predominantly couple to the 
vector interpolator, and the other ones would couple pre- 
dominantly to the pseudotensor operator. Asymptoti- 
cally each of these mesons would belong entirely to one 
of the two representations, (0, 1) -I- (1,0) or (1/2, 1/2)^, 



and could not be created by the operator that transforms 
according to the other representation. 

In lattice simulations, using the variational method, it 
is possible not only to reliably separate a given excited 
state and measure its mass [ISj] but also to define and 
measure a ratio of couplingsof different lattice operators 
to the given excited state [l9|. The two chiral represen- 
tations (0, 1) -I- (1, 0) and (1/2, 1/2)^ form a complete and 
orthogonal basis (with respect to the chiral group) for the 
p-meson. Consequently using the variational method we 
are able to study a mixing of two representations in the 
excited p-meson and see whether or not a given excited 
p-meson at low resolution (infrared) scale couples pre- 
dominantly to one of the representations and decouples 
from the other. 

Assuming that the set of interpolating operators Oi (t) 
is projected to vanishing spatial momentum, the energies 
of the states S*^"-' and the coefficients giving the overlap 
of operators with the physical state, a^" — {0\Oi\n) , can 
be extracted from the cross-correlation matrix 



C(t),, = (O.(t)Ot(0))=^af)a(")*e 



(n) (n)* -E^"h 



(1) 



It can be shown [l8| (for a recent discussion see 26l.l27|) 
that the generalized eigenvalue problem 



a(i),,«^"' = A(")(i,<o)5(io).j^i$"^ 



(2) 



allows to recover the correct eigenvalues and eigenvec- 
tors within some approximation. At t = to all eigenval- 
ues are 1 and the eigenvectors are arbitrary. Energies of 
the subsequent states can be extracted from the leading 
exponential decays of each eigenvalue. 

Ratios of couplings of the different operators to the 
physical states can be obtained as (cf., [l9|) 



C{t)ijU 
C{t)kj 
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(3) 



A decomposition of a hadron depends on the resolution 
scale, i.e., what we see in our microscope, depends on 
its resolution. If one uses point-like lattice interpolators 
then the resolution scale 1/i? is determined by the lattice 
spacing a. We are interested to study a decomposition 
of a hadron at a very low resolution scale, determined by 
the hadron size. For that we cannot use a large a because 
a proper matching with the ultraviolet (continuum) be- 
havior of QCD will be lost. Given a fixed, reasonably 
small, value for a, a small resolution scale 1/ R can be 
achieved by the g aug e-invariant smearing of the point- 
like interpolators [19j. We use the interpolator smeared 
over the size R in physical units such that R/a ^ 1, so 
even in the continuum limit a — >■ we probe the hadron 
structure at the scale fixed by R. Changing the smear- 
ing size i?, we can study the hadron content at different 
resolution scales of the continuum theory at a — >■ 0. For 
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FIG. 1: The vector meson mass mv is plotted against m„ for 
all three sets. Black circles represent the ground state, p, and 
red squares represent the first excitation, p' . The experimen- 
tal values are depicted as magenta crosses with decay width 
indicated (color online). 



TABLE I: Specification of the data used here; for the gauge 
coupling only the leading value Plw is given, mo denotes the 
bare mass parameter of the CI action. Further details on the 
action, the simulation and the determination of the lattice 
spacing and the n- and p-masses are found in 24, ^] . 

Set Plw a mo #conf a [fm] m^ [MeV] nip [MeV] 



A 4.70 -0.050 
B 4.65 -0.060 
C 4.58 -0.077 



200 
300 
300 



0.1507(17) 
0.1500(12) 
0.1440(12) 



526(7) 911(11) 
469(5) 870(10) 
323(5) 795(15) 



this purpose we use Jacobi smearing (281 ) , which provides 



a Gaussian shape of interpolators in spatial directions 
of size R at both source and sink. Such definition of 
the resolution scale is similar to the experimental defini- 
tion where a resolution is determined by the momentum 
transfer in spatial directions. 

3. Lattice simulation and results. In our excited 
hadron s pec troscopy program, first for quenched config- 
22L I29I [SOJ and recently for dynamical fermions 



urations 



24 Isij the Liischer- Weisz gauge action is used [32] • For 
the fermions the Chirally Improved (CI) Dirac operator 
is adopted, which has better chiral properties than the 
Wilson Dirac operator |33|. 

In this study of the chiral decomposition of the excited 
p-meson we use three sets of dynamical configurations 
with two mass-degenerate light sea quarks on lattices of 
size 16^ X 32 (see Tab.©. 

The following set of four interpolators is used, 

Oi = Unl'^dn , O2 = UwYdw , (4) 

O3 = Un"f*Y<in , Oa = Uwj'^^^dio . (5) 

7* is one of the spatial Dirac matrices, 7* is the 7-matrix 
in (Euclidean) time direction, and the subscripts n and w 
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FIG. 2: The ratio av/ar is plotted against the smearing 
width R for all three data sets. Black circles represent the 
ground state and red squares the first excitation. Broken 
lines are drawn only to guide the eye (color online). 



(for narrow and wide) denote the two smearing widths, 
R « 0.34 fm and 0.67 fm, respectively. This bunch of 
operators allows one to extract both the ground and the 
first excited state of the p-meson, see Fig. [T] 

In Fig. [2] we show the i?-dependence of the ratio ay far 
both for the ground state p-meson and its first excited 
state. This ratio of the vector to pseudotensor coupling 
gives us the information on the chiral decomposition of 
both states in terms of the representations (0, 1) -|- (1, 0) 
and (1/2, 1/2);,. We observe that this ratio for the ground 
state at the smallest possible resolution scale (largest 
smearing radius) R — 0.67 fm approaches a value close 
to 1.2. This implies a strong mixture of both represen- 
tations in the p-meson wave function. Using a unitary 
transformation from the chiral basis to the ^^^^Lj basis 



^^1) = 



'D,) = 



-|(0,1) + (1,0);1 — ) + 
^ 1(1/2, l/2)fc; r 

^|(0,1) + (1,0);1 — )- 
^|(l/2,l/2)fc;r 



(6) 



one obtains that the ground state p is a '^ 5*1 state with a 
tiny admixture of the ^Di wave, O.ggTl^S"!) -0.073pDi). 
This implies that at a resolution fixed by the p-size the 
p(770) is approximately a '^S'l state in agreement with 
the quark model. 

Our main result in the present report is that the chi- 
ral decomposition of the first excited state, p' — p(1450), 
and its scale dependence is very different. In this case at 
large R a leading contribution comes from the (1/2, l/2)h 



representation. Given only two different R values we can- 
not reliably extrapolate to the scale of 1 fm, as suggested 
by the size of excited rho-meson. It is clearly seen, how- 
ever, that such a ratio at i? '^ 1 fm is very small; it can 
be both positive and negative in sign. 

This behavior has a clear interpretation. In the deep 
ultraviolet one expects from the conformal symmetry of 
QCD that the pseudotensor interpolator decouples from 
the /5-mesons. Consequently towards small R the ratio 
av/cLT niust increase. At large R the ratio determines 
a degree of the chiral symmetry breaking in the infrared 
region, where mass is generated. One observes that such 
a breaking for the p(1450) is insignificant, in contrast 
to the p(770). Since the chiral decomposition of the p' 
state is dominated by one of the chiral representations, 
it indicates a smooth onset of effective chiral restoration. 
Given that this leading representation of p' is (1/2, l/2)h 
one predicts that in the same energy region there must 
exist a hi (and not a oi) meson. Inspecting the Particle 
Data Group [j| one finds that there is indeed the state 
/ii(1380) and no ai in the same energy region. 

One naturally arrives at the following identification 
for /5-mesons. Chiral symmetry is strongly broken in 
the p(770), since it is a strong mixture of (0, 1) + (1,0) 
and (1/2, l/2)(,, and consequently its "would-be chi- 
ral partners" have a much larger mass: ai(1260) and 
/ii(1170). AU these lowest-lying states p(770), ai(1260) 
and /ii(1170) cannot be assigned to any chiral represen- 
tation. In the first excited state, /9(1450), a contribu- 
tion of (1/2, l/2)f, is much bigger than a contribution of 
the other representation, and consequently its approxi- 
mate chiral partner is the /ii(1380). Then, the next ex- 
cited p-meson, p(1700), is most probably dominated by 
(0, 1) -|- (1, 0), which is supported by existence of a nearby 
ai state, ai(1640). Note that there is no room for this 
ai(1640) meson within the quark model [3|, |j|. 

While we do not know the precise value of the ratio 
av/dT for the p(1450) at i? ~ 1 fm, it is indicative that 
this ratio is very small. Then we can qualitatively esti- 
mate the angular momentum content of p(1450) in the 
infrared. Assuming a vanishing ratio, one obtains the 
following partial wave content, ^/T/3\^Sl) — ^/2/3\^Dl). 
A significant contribution of the '^Di wave is obvious. A 
possible variation of the ratio at large R changes slightly 
numbers for the partial wave decomposition, but does 
not change the qualitative result. This result is not con- 
sistent with p' to be a radial excitation of the ground 
state p-meson, i.e., a ^^i state, as predicted by the quark 
model. 
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